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Abstract. We propose a segment description for turbulent solutions to evo- 
lution equations in an effort to develop an effective description rather than the 
classical Reynolds average. The new description has connections with symbolic 
dynamics and shadowing technique. The challenge of future study is how to 
effectively implement the description numerically. 



1. Introduction 

Chaos and turbulence have no good averages [3]. The matter is more funda- 
mental than just poor understanding of averages. The very mechanism of chaos 
leads to the impossibility of a good average . From what we learn about chaos 
in partial differential equations [2], turbulent solutions not only have sensitive de- 
pendences on initial conditions, but also are densely packed inside a domain in the 
phase space. They are far away from the feature of fluctuations around a mean. 
In fact, they wander around in a fat domain rather than a thin domain in the 
phase space. Therefore, averaging makes no sense at all. One has to seek other 
descriptions. In this article, we present a description - named segment description. 
Segment description approximates a turbulent solution by an infinite sequence of 
segments over an infinite time interval. The total number of different segments is 
finite, and all the segments are over a fixed finite time interval. Thus, each segment 
can be calculated rather accurately. The segment description also has connections 
with both symbolic dynamics and shadowing technique [21 . But, our emphasis here 
is upon applications rather than mathematical abstract arguments. In our view, 
there are two most significant problems in the study of turbulence: 1. an effective 
description of turbulent solutions rather than the Reynolds average, 2. control of 
turbulence. Here the segment description is aiming at problem 1 by offering an 
alternative description. 

In a sense, the Reynolds average is based upon the feature of fluctuations 
around a mean, which turbulent solutions do not have. Consider the Navier-Stokes 
equations in the form: 

(1.1) Ui^t + UjUi^j = -p^i + Rc""^ Ui^jj + fi , Ui^i = ; 
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where Ui's are the velocity components, p is the pressure, fi's are the external force 
components, and Re is the Reynolds number. According to the Reynolds average, 
one splits the velocity and pressure of fluids into two parts: 

Ui ^ Ui + Ui , p = P +p 

where the capital letters represent relatively long time averages which are still a 
function of time and space, and the tilde- variables represent mean zero fluctuations, 

- {u,) , (w,) = , P = (p) , (p) = . 

A better interpretation is by using ensemble average of repeated experiments. One 
can derive the Reynolds equations for the averages, 

(1.2) Uu + U-iU^^j = -P. + RC-I C/.jj - {U^UJ)^J + h , U^^^ = . 

The term {uiUj) is completely unknown. Fluid engineers call it Reynolds stress. 
The Reynolds model is given by 

(1.3) {u^u,) = -R-^ U,,j , 

where i? is a constant. There are many more models on the term (uiUj) [1]. But 
no one leads to a satisfactory result. One can re- interpret the Reynolds equations 
(|1.2p as control equations of the orginal Navier-Stokes equations ()1.1|) . with the 
term {uiUj)j being the control of taming a turbulent solution to a laminar solution 
(hopefully nearby) [3]. The Reynolds model (|1.3|) amounts to changing the fluid 
viscosity which can bring a turbulent flow to a laminar flow. This laminar flow 
may not be anywhere near the turbulent flow though. Thus, the Reynolds model 
may not produce satisfactory result in comparison with the experiments. Fluid 
engineers gradually gave up all these Reynolds' type models and started directly 
computing the original Navier-Stokes equations (jl.ip . 

2. Segment Description 

Let S be the phase space which is a metric space, on which a flow P* {t > 0) is 
defined. F* is C° in t, and for any fixed t E (0, +oo), F* is a map. For instance, 
S can be a Sobolev space defined on the 2-tori, and F* is the 2D Navier-Stokes flow. 
Let A he a compact absorbing set on which the dynamics of F* can be turbulent. 
The attractor is a subset of A. A can be chosen to be just the attractor. But a 
compact absorbing set is easier to obtain than an attractor in numerics. Let B be 
a ball that includes A, such that 

(2.1) \\DF\x)\\<C{t), yxeB, 

where C{t) is a continuous function of i G [0, -l-oo). Let D be any subset of S, the 
diameter of D is defined as 

diamcterjl?} ~ sup distance(a::, y) . 

x,yeD 

Let Bs{x) be the ball centered at x with radius 5 > Q. Fix a F > and a small 
e > 0, for any x E A, there exists a (5 > [e.g. 6 = ce for some c > by (|2.ip ] such 
that 



(2.2) 



sup diameter {F* (5^(2;))} < e 
te[o,T] 
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The union of all such balls 

U Bsix) 

is a covering of A. Since A is compact, there is a finite covering 

U Bs„{Xn) ■ 
l<n<N 

We can choose this covering to be minimal. From this covering, we can deduce a 
mutually disjoint covering: 

BsSxn) = S5„(x„)f| y Bs^{x„,) , (n = I,--- ,Ar-l) ; 

n+l<m<N 

Bs^{xn) = Bs^{xn) ■ 
Then lJi<ri<A' Bs„{xn) is a mutually disjoint covering of A. Thus 

^Ar, ^ Bsjxn) C\ A , (l<n<7V)} 

form a partition of A. Let /i be a Borel measure defined on A, then the metric 
entropy of the partition is 

N 

ij = - ^/x(A„)iog^(AO • 

In this article, we are not going to emphasize either metric entropy or topological 
entropy [4] since they are essentially sums of positive Liapunov exponents and Lia- 
punov exponents are very easy to compute in numerics. Unlike the usual partitions, 
here A„'s are of small size. On the other hand, we do not take the size to the zero 
limit [1] in which case we will be led to Liapunov exponents. Our goal is to obtain 
a uniform approximation of the turbulent solutions over an infinite time interval. 

Definition 2.1 (Segments). Define the segments s„ as follows 

s„ = U F\xn) , (n = I,-- - ,iV) . 

t6[0,T] 

Definition 2.2 (Segment Description). For any x e ^, the orbit 

(2.3) ^(t)^{F\x), ie[0,+oo)} 
has the segment representation 

(2.4) 77(i) = {s„oS„,s„, •••} , te[0,+cx)) 
where 

F'^{x)(,An,, Tij e {I,--- ,iV} , 0<j<+w. 
We call ri{t) an e-pseudo-orbit. 

Lemma 2.3. For any t e [0, +oo), 

distance {^{t),r]{t)) < e . 

Proof. The lemma follows from the condition (|2.2p and the Definition 12. 21 □ 
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We define a maximal difference function as follows 

Md{t) — max distance (s„i(i), s„(t)) . 

l<m,n<N 

If the absorbing set A is actually an attractor, then Md{t) should be almost constant 
in t and almost the diameter of the attractor. In general, Md(t) serves as a measure 
of the shrinking of the diameter of the absorbing set A under the flow. 

2.1. Markov Description . For any 1 < m, n < A^, if 

(2.5) F^{A^)nA„^IJ} , 

then we assign the value r,„„ — 1; otherwise F^n = 0. These F„i„'s form a N x N 



matrix F called the transition matrix. In the segment representation (12. 4|) . the first 
segment s„q can be any of the N choices, while the second segment s„j is admissible 
only when F„„„^ = 1. The admissibility of the third segment s„2 depends on both 
Sno a-nd Sni- In general, the admissibility of the j-th segment Snj_i depends on 
all segments in front of it s„q ■ • ■ s„ -_.j. Thus we are dealing with a non-Markov 
process. Nevertheless, we believe that the Markov description is still significant. 
The Markov description extended the admissibility by the criterion that for any 
j = 1, 2, • • • ; is admissible if F„^_j„^. — 1. We can define a probability as follows 

(2.6) Pmn = M (F^iAm) O A„) ^ fl{F^ {Am)) , 

where fi is the Borel measure defined on A. 

Definition 2.4. Let Si be the space consisting of elements of the form 

a = {nonin2 ■ ■■} 

where no S {!,•■• , N} and the admissibility of nj is decided by F„j_j„j = 1 for 
j = 1, 2, • • ■ . The product topology in Si is given by taking as the neighborhood 
basis of any element 

* f * ^ ^ 1 

a |no71i7l2 • • • } 

the cylinder sets 

Wj = {a e Si I rik = nl {k < j)} 

for j — 1, 2, ■ • • . The left shift Xi is defined by 

Xi{a) = {nin2n3 • ■ • } , Va = {no?T-in2 • • • } ■ 

The pair (Si, xi) is called a subshift of finite type. 

Lemma 2.5. // each row ofT has more than one entry to be 1, then the left 
shift xi has sensitive dependence on the initial conditions. 

Proof. For any large j, let a, a £ Wj and 

a = {fiQnin2 • • ■ } , a = {no"i»^2 • • • } ; 

then 

fik = nk {k < j) ■ 

Since each row of F has more than one entry to be 1, we can choose such a and a 
that fij ^ fij. After j iterations of xi, we have 

X{{a) = {njnj+inj+2 • • ■ } , xi(a) = {^j^i+i"j+2 • ■ • } ■ 
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Thus 

xi(a)^T4^i ofxi(5) . 

That is, after j iterations of xi, the two elements are far away. The proof is 
complete. □ 

In the case that not all the rows of F have more than one entry to be 1, the 
left shift xi may exhibit transient chaos (intermittence) or regular dynamics. 

Definition 2.6. Let S2 be the space consisting of elements of the form 

where no G {1, ■ ' ' 1-^} ^-nd the admissibility of rij is decided by r„j_j„^. = 1 for 
j = 1,2, • ■ • . The product topology in Si induces a product topology in S2 with 
the neighborhood basis of any clement 

the cylinder sets 

V, = {((t) eE2\nk = nl (fc < j)} 
for j = 1, 2, ■ • • . The left shift X2 is defined by 

X2(C(0) = {SniS„2S„3 • ■ • } , VC(t) = {SnoSmSn^ ■■■} . 



Definition 2.7. Let S be the space consisting of orbits starting from any point 
in A 

^{t)^{F\x), te[0,+oo)}, VxeA. 

The product topology in S2 induces a topology in S via the map (p : E S2 
which maps (|2.3p to (|2.4p . For any ^*{t) € S, ^*{t) has a segment representation 
as given in (|2.4p 

7y*(t) = |s„.s„-s„* • ■ • } , iG[0, +00) 

at which the neighborhood basis V^ 's are defined in Definition 12.61 The preimages 
of the Vj 's under Lp form the induced neighborhood basis at ^* (t) 

Uj^V^-HVj) (j = l,2,-..). 



In general, the map (p : E S2 which maps (|2.3p to (|2.4p is neither 1-1 nor 
onto. Under the topologies defined in Definitions 12.61 and 12. 7[ p is a. continuous 
map. 

Theorem 2.8. There is a continuous map (j) : E k-s- Si such that the following 
diagram commutes: 

E ^ Si 
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Proof. Let ip be the homeomorphism 

-0 : S2 1-^ Si , ^p{{snoSnlS„^ ■■■}) = {nonin2- ■■} . 

Then 

(j) = ip o (p : Si— >Si 

is a continuous map. For any ^(t) G S, ^(t) has a segment representation given by 
Then F'^{^{t)) has the segment representation 

Thus the diagram commutes. □ 

Theorem 12 .81 serves as a Markov description of (S, F^). Notice that the proba- 
bihty (|2.6p generates a Markov measure on the subshift of finite type (Si, xi)- The 
Markov measure of the cyhnder sets Wj in Definition 12.41 is given by 

The Kolmogorov-Sinai entropy is 

N N 
m— 1 n— 1 

Using the Markov description, one can estimate interesting quantities. Let Q{x) 
be an interesting quantity (for example dissipation rate or energy for Navier-Stokes 
equations), one can calculate 

< = sup Qisnit)) , ql = inf Q(s7.(t)) , (n = 1, • • • , TV) . 

te[o,T] te[o,T] 

For any finite time interval [0,mr] (m e Z+), let G ^7iq, then one can track all 
the possible finite sequences starting from uq 

(2.7) {"■0"'l"2 • • -"m-i} 

by the condition r7i =1 (j = 1, • • • , m — 1) (cf. Definition 12. 4p . fn particular, 
one can identify the subset I™ of {1, ■ • ■ , N} that consists of elements appearing 
in the sequences (I2.7p . Then one can estimate Q{£,{t)) {£,{0} = xq) over the time 
interval [0,mT] by 

= max q!^ , = max . 

That is 

q'<Q{m)<q\ ie[o,mr]. 

One should expect that numerically it should be quite easy to track the sequences 
(|2.7p up to very large m. 
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2.2. Non-Markov Description. In general, not all the elements in S2 can 
be realized by an orbit in E. Thus all the orbits in S are labeled by a subset of S2 
via (|2.4p . To better describe the segment description (|2.3p - (|2.4p . we introduce the 
transition tensors 

r r 

where Tngm is defined by p.Sp . if 

then r„Q„j„2 = 1, otherwise Tnonin-i — 0. Similarly one can define the rest of the 
transition tensors. 

Definition 2.9. We say that the transition tensors are expanding if for any 
^noni---nj — 1, there is an integer m > such that there are at least two different 
values of n^+rn for which r„„„j...„^....„^^^ = 1. 

Theorem 2.10. // the transition tensors are expanding, then : E i— > E 
has sensitive dependence on initial conditions. That is, F'^ is chaotic. 

Proof. For any ^(t) e S and any neighborhood Uj fDefinition l2.7p . since the 
transition tensors are expanding, there is another element ^* (t) G Uj such that the 
segment representations of ^ (t) and ^* (t) are of the form 

{ -^riQ "^ni '^712 ' * ' "^Tij ' ' ' •^rij^Tn ' ' ' } and Stjq Syj^ ^7^2 * ' * -^Tij ' ' ' ' * ' ^ 

where nj+m ^ n*^„^. Then 

F<-J+"^^T{^*{t)) ^ Ui of F(J+™)^(C(t)) . 
Thus F'^ has sensitive dependence on initial conditions. □ 

If the transition tensors are not expanding, then F"^ may exhibit transient 
chaos (intermittence) or regular dynamics. 

Definition 2.11. Let be the space consisting of elements of the form 

a = {noni?i2 • • • } 

where for any j = 0, 1, 2, ■ • ■ , G {1, • ■ • , N}; and r„Q...„^ = 1 for any fc = 1, 2, • • • . 
The product topology in and the left shift x are defined as in Definition 12.41 

It is easy to see that is the image of cj) in Theorem 12.81 Thus as a corollary 
of Theorem 12. 8[ we have the following theorem. 

Theorem 2.12. There is a continuous map cf) : E i— > which is onto, such 
that the following diagram commutes: 
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3. Numerical Implementation 

In a numerical implementation, one starts with a compact absorbing set A (for 
example, a ball in a subspace spanned by a finite number of Fourier modes) in a 
Banach space. Then one can collocate as many points as possible in the absorbing 
set A. Starting from these points, one can generate the segments over a time 
interval 

s„(t) , (n = I,--- ,iV) , te[o,r]. 

One can also compute the maximal difference function 

Mdit)= max \\Smit) - Sn{t)\\ . 
l<m,n<N 

To implement a Markov description, one can compute the norm of the gradient of 
restricted to A at s„(T): 

Pn = ||VF^(s„(r))U|| . 

Let 

rn 

and choose such that 

\\sn{T) - s„.(0)|| - min {||s„(T) - Sfc(0)||} . 

Then any Sm(0) (1 < to < N) inside the ball Bp^r„ (sn, (0)) implies that Sm is 
admissible to follow Sn in the segment sequence (cf. Definition 12. 6p . One can also 
calculate interesting quantities as described at the end of subsection 12. II 

The challenge in the numerical implementation of the segment description lies 
at the dimension of the absorbing set A. Along each dimension, if we collocate 10 
points, and if the dimension of ^ is n, then we have total 10" collocation points. So 
when the dimension n is greater than 9, the computation will be very challenging. 
The hope is that future computer technology may overcome this difficulty. 

For low dimensional system like Lorenz equations, the numerical implemen- 
tation of the segment description should not be difficult. For 2D Navier-Stokes 
equations with periodic boundary condition, numerical implementation of the seg- 
ment description for a Galerkin truncation of less than 6 modes should not be 
difficult either. 

4. Conclusion 

In this article, we introduced a segment description of turbulent solutions. The 
goal we are aiming at is to develop an effective description of turbulent solutions 
rather than the Reynolds average. Via the segment description, every turbulent 
solution can be approximated uniformly over an infinite time interval by a pseudo- 
solution built by segments of solutions. In dynamical systems, shadowing lemma 
addresses the converse of the above claim. In certain settings, shadowing lemma 
holds [2] and for every pseudo-solution, one can find a true solution nearby. Even 
though the total number of segments in the segment description is finite, the number 
can be very large when the dimension of the phase space is large. This poses a 
challenge in numerical implementation of the segment description. 



in {P„(0)-s,(0)||} , 
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